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Abstract 

To understand the nature of supercritical accretion, we systematically analyze the RXTE/PCA data 
of GRS 1915+105 in its quasi-steady states, by choosing data with small variability during 1999 - 2000. 
We apply a multicolor disk plus a thermal Comptonization model and take into consideration accurate 
interstellar absorption, a reflection component (with an iron-K emission line), and absorption features from 
the disk wind self-consistently. The total luminosity ranges from ~ 0.2Le to slightly above Le- There is a 
strong correlation between the inner disk temperature and the fraction of the disk component. Most of the 
Comptonization-dominated (> 50% total flux) spectra show T- ln ~ 1 keV with a high electron temperature 
of > 10 keV, which may correspond to the very high state in canonical black hole X-ray binaries (BHBs). 
By contrast, the disk-dominated spectra have T ln ~ 2 keV with a low temperature (< 10 keV) and optically 
thick Comptonization, and show two separate branches in the luminosity vs. innermost temperature (L 
- Ti n ) diagram. The lower branch clearly follows the L oc X^-track. Furthermore, applying the extended 
disk blackbody (or p-free disk) model, we find that 9 out of 12 datasets with disk luminosity above 0.3Le 
prefer a flatter temperature gradient than that in the standard disk (p < 0.7). We interpret that, in the 
lower branch, the disk extends down to the innermost stable circular orbit, and the source is most probably 
in the slim disk state. A rapidly spinning black hole can explain both the lack of the L cx T^-track and 
a high value of spectral hardening factor (~ 4) that would be required for a non-rotating black hole. 
The spectra in the upper branch are consistent with the picture of a truncated disk with low temperature 
Comptonization. This state is uniquely observed from GRS 1915+105 among BHBs, which may be present 
at near-Eddington luminosity. 
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1. Introduction 

The dynamics and emission properties of supercriti- 
cal (or super-Eddington) accretion flows have been long- 
standing issues since the 1970s (Shakura and Sunyaev 
1973). Theoretical studies on supercritical accretion are 
far from being complete due to the complex interaction 
between radiation and matter, despite recent progress in 
numerical simulations (e.g. Ohsuga et al. 2009). From 
observational studies, supercritical accretion appears to 
be present in some extragalactic objects, such as Narrow- 
line Seyfert 1 galaxies (NLSls) and Ultraluminous X-ray 
sources (ULXs) (Kato ct al. 2008, and references therein). 
In some NLSls, the black hole mass can be estimated 
by means of reverberation mapping (e.g. Wandel et al. 
1999; Kaspi et al. 2000). In the case of ULXs, however, 
the black hole mass is unknown and thus remains contro- 
versial. Nevertheless, their luminosity is higher than the 
Eddington luminosity of a stellar mass black hole. Black 
hole mass is the key to initially identifying a supercritical 
accreting system since the Eddington luminosity depends 
only on the black hole mass. In order to do a comprehen- 
sive study on supercritical accretion processes, we need an 
object whose black hole mass can be well-constrained. In 
addition, we also need adequate observational data. 



One of the best objects for the study of supercritical 
accretion is the Galactic microquasar, GRS 1915+105. 
GRS 1915+105 has been active since its discovery by 
the WATCH all-sky X-ray monitor on GRANAT (Castro- 
Tirado ct al. 1992; 1994) on 1992 August 15. Its distance 
and inclination angle have been estimated from its super- 
luminal jets (Mirabel and Rodriguez 1994) to be 12.5 + 1.5 
kpc (see also Chaty et al. 1996; Dhawan et al. 2000) and 
70° ±2°, respectively. Dozens of types of peculiar vari- 
ability in the X-rays have been observed (see Belloni et al. 
2000; Klein- Wolt et al. 2002; Hannikainen et al. 2005), 
most of which have never been observed in other Galactic 
X-ray sources, including other known microquasars. GRS 
1915+105 is thus known to be a peculiar X-ray source 
in our Galaxy. The jet feature and the variability of GRS 
1915+105 are believed to be an important key to studying 
the coupling between jets and accretion disks in general 
(e.g. Mirabel et al. 1998; Neilsen and Lee 2009). 

While optical observations are severely hampered since 
GRS 1915+105 is located close to the Galactic plane at 
I = 45.37, b = -0.22 (Castro-Tirado et al. 1994), observa- 
tions at other wavelengths prove to be quite fruitful. From 
infrared spectroscopy, GRS 1915+105 was suggested to be 
a low-mass X-ray binary (Castro-Tirado et al. 1996). This 
result is further supported by Greiner et al. (2001a) and 
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they also suggested that the donor star is a K - M giant 
with a narrow mass range ~ 1.0 - 1.5 AT©. The dynamical 
mass estimation results in a compact object of 14 ± 4M 
(Greiner et al. 2001b). This classifies GRS 1915+105 as 
the most massive black hole X-ray binaries (BHBs) in our 
Galaxy (for extragalactic BHBs, see, e.g. Orosz et al. 
2007; Prestwich et al. 2007). 

Since the dynamical mass can be estimated, GRS 
1915+105 is known to prefer being in a state where the lu- 
minosity is close to or even exceeding its Eddington lumi- 
nosity (Done et al. 2004), which corresponds to the super- 
critical accretion process. Provided that GRS 1915+105 
docs really exhibit supercritical accretion, it can be used 
as a guide for the interpretation of other supercritical ac- 
creting objects that lack sufficient data for a comprehen- 
sive study. 

We use archival RXTE data mainly from 1999 to 2000 
and apply the so-called disk-corona model to the extracted 
spectra. Our focus on the study of the supercritical ac- 
cretion in GRS 1915+105 distinguishes this work from 
other previous extensive studies (e.g. Done et al. 2004; 
McClintock et al. 2006; Middleton et al. 2006). We also 
need to stress that we do not attempt to explain the 
nature of GRS 1915+105 with its richness of variabil- 
ity on timescales from years down to milliseconds (see, 
e.g. Morgan et al. 1997; Belloni et al. 2000) and its non- 
canonical soft and hard states. Those studies have been 
done by many authors (see, e.g. Belloni et al. 1997a, b) 
and the most recent work can be found in Rodriguez et al. 
(2008a, b). On the contrary, we choose data from when 
the source is more or less stable, without strong variabil- 
ity features as classified in Belloni et al. (2000). The 
constraint on the mass of GRS 1915+105 from dynamical 
estimates helps us to study the detailed relation between 
the X-ray spectra and mass accretion rate in terms of its 
luminosity during its supercritical accretion state. 

The layout of this paper is as follows: We describe our 
data sample in section 2, and continue with describing the 
model and the fitting analysis in section 3 and 4, respec- 
tively. Section 5 is devoted to the discussion, and finally 
section 6 concludes the paper. Throughout our study we 
adopt the following values: a distance of 12 kpc, an incli- 
nation angle of 70°, and a black hole mass of 14M Q . 

2. Data 

For the purpose of the present study, we deliberately 
choose the data when the source was bright and more or 
less stable, that is when no large flux variations are ob- 
served. Therefore, we only choose data from classes 0, 8, 
and x from the 12 known classes of X-ray variability in 
GRS 1915+105 (Belloni et al. 2000). The variability selec- 
tion is made by eye from the lightcurve available through 
the RXTE standard products. 

GRS 1915+105 has been monitored by RXTE since 
1996, once a week on average. We select the RXTE data 
of GRS 1915+105 between January 1999 and May 2000. 
This range of observation times contains data from epochs 
3B, 4 and 5, and consists of a total of 216 datasets to be 
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Fig. 1. RXTE/ ASM lightcurve. The minimum count for soft 
state detection, 30 cts/s, is marked with the dotted line. The 
arrows mark the datasets of our final sample (see text). 

examined. Focusing on data with small variability, we 
manage to keep 23 datasets for further fitting analysis. 
In addition, we also analyze two disk-dominated datasets 
from 1997 (epoch 3A) and 1998 (epoch 3B), previously 
studied by Done et al. (2004). To sum up, our final sam- 
ple consists of 25 datasets (epoch 3A, 3B, 4, and 5). The 
details of the selected datasets are given in Table 1. 

We follow the standard data reduction procedure for 
RXTE data. We use PCA data only for the fitting en- 
ergy range we are interested in, 3-25 keV, with 16-s 
time-resolution of the Standard 2 data. The lightcurves 
and the spectra of the source and background were ex- 
tracted using SAEXTRCT version 4.2e. The latest version 
of RUNPCABACKEST (version 4.0) was used to estimate the 
background. The response files were created by using 
PCARSP version 10.1. The good time intervals are selected 
using the following conditions: elevation angle greater 
than 10 degrees, the OFFSET less than 0.02, and we use 
all detectors which were in operation during the obser- 
vation time. Dead-time correction is applied and 1.5% 
systematic error is included in each spectral bin. We find 
no significant change in the fitting parameters when a 1% 
systematic error is applied. Figure 1 shows the one-day 
average RXTE/ASM lightcurve of GRS 1915+105 (sum 
band intensity) from MJD 50700 - 51700. The arrows 
mark the datasets of our final sample. 

We follow Belloni et al. (2000) for the description of 
states A, B, and C in Table 1. States A and B correspond 
to disk-dominated spectra while state C corresponds to 
the power-law type spectra. State A has lower flux and 
disk temperature (~ 1 keV) compared to those of state B. 
State B is associated with high mass accretion rate and 
its spectra resemble that of the very high state (VHS) in 
normal BHBs (e.g. Belloni et al. 2000, and the references 
therein). State C is mostly found in variability class \i 
while states A and B can be found in variability class 8. 
Variability class cf> mostly consists of state A. 



GRS 1915+105 at Near-Eddington Luminosity 

Table 1. The final sample. 



3 



Dataset 


/~\1, - „ i • _ 1 1 \ 

(Jbscrvation ID 


(Jbscrvation Date 


Observation lime 


Duration (s) 


Glass (otatej 


1 


on i no rr A A 

20402-01-55-00 


1 A A T 11 IT 

1997-11-17 


A ET AT A/ 1 

05:27:06.6 


15378 


5 (B,(A)) 


2 


30703-01-08-00 


1 a a a a 1 a 

1998-02-14 


AO P" 1 r" p' 

23:51:55.3 


4868 


1 / A \ 

(A) 


3 


atao ai 1 /i ai 

40703-01-14-01 


1 CI ("1 CI A r* AO 

1999-05-08 


/~1 A O IO T 

04:33:18.7 


4319 


c /t-> /AW 

5 (B,(A)) 


4 


a a t a o ai 1 ^ a o 

40703-01-14-02 


1 AAA A r* AO 

1999-05-08 


A/"* AA 1 A 

06:09:19 


3085 


e /-r> / A \ \ 

5 (B,(A)) 


5 


^ A17A0 ai a a a a 

40703-01-20-02 


1 AAA AT AO 

1999-07-03 


A 1 P* 1 O T T 

01:51:37.7 


4190 


X (C) 


6 


a a t a o ai 00 ai 

40703-01-22-01 


1 AAA AT O "I 

1999-07-21 


A 1 P* T A A P" 

01:57:09.5 


2699 


X (C) 


7 


^ ATAO r\ 1 AA 

40703-01-28-00 


1 AAA AO AO 

1999-08-28 


1 A OA AA 

18:30:00.6 


2891 


5 (B,(A)) 


8 


a ata a 1 00 a 1 

40703-01-28-01 


1 AAA AO AO 

1999-08-28 


A A 1 T 1 1 O 

20:17:11.3 


2220 


c /-r> /AW 

<5 (B,(A)) 


9 


h /"a t a ai oa a a 

40703-01-29-00 


1 AAA AA A f" 

1999-09-05 


IO 1 A A A A 

18:12:41.2 


3210 


e /t-> /AW 

<5 (B,(A)) 


10 


a n^no ai nn ai 

40703-01-29-01 


1 AAA AA A f" 

1999-09-05 


1 A r A 1 1 

19:59:11.3 


2520 


e /t-> /AW 

«5 (B,(A)) 


11 


<ni 1 c ai no n 1 

40115-01-08-01 


1 AAA 1 A 1 A 

1999-10-19 


11 r* r* A /"* 

11:55:03.6 


3514 


r /t-> /AW 

<5 (B,(A)) 


12 


yi atao a 1 o^i a a 

40703-01-34-00 


1 AAA 1 A AA 

1999-10-20 


1 FT 1 A A C P* 

15:19:06.5 


2731 


e / ~r~) /AW 

<5 (B,(A)) 


13 


< A T A O ai O/i a 1 

40703-01-34-01 


1 AAA 1 A AA 

1999-10-20 


1 T A A or - 1 

17:09:35.1 


1862 


£ /-r> /AW 

<5 (B,(A)) 


14 


( ATAO r\ 1 or" AA 

40703-01-35-00 


1 AAA 1 A C\f 

1999-10-25 


IP" AT AO P* 

15:07:08.5 


3108 


X (C) 


15 


/I ATAO A 1 O /"* A A 

40703-01-36-00 


1 AAA 1 1 AA 

1999-11-02 


1/1 /i O O T A 

14:48:37.4 


2415 


X (C) 


16 


a a^ao a 1 OO A 

40703-01-38-03 


1 AAA 11 1 r 

1999-11-15 


1 A AA AO O 

14:09:02.3 


3420 


X (C) 


17 


40703-01-38-02 


1 AAA 11 1 

1999-11-15 


IP" A O A 1 ^ 

15:43:01.6 


3540 


X (C) 


1 s 


in7n"? m ni 

4U / Uo-Ul-oo-Ul 


1 QQQ 11 1 ^ 

1-10 


17.1q.n1 r 

1 / .lo.Ul.D 


oyuu 


X (M 


19 


40703-01-39-00 


1999-11-21 


12:09:00.9 


3901 


S (B,(A)) 


20 


40703-01-40-00 


1999-12-01 


11:45:03.2 


3900 


X(C) 


21 


40703-01-40-02 


1999-12-01 


14:56:22.3 


2912 


X(C) 


22 


40703-01-40-03 


1999-12-01 


16:30:03.5 


3240 


X(C) 


23 


40702-01-04-00 


2000-02-05 


05:58:34.5 


3004 


x(c) 


24 


40403-01-10-00 


2000-02-14 


10:28:11.1 


2677 


x(C) 


25 


50703-01-10-01 


2000-05-17 


03:15:10.7 


2753 


X(C) 



3. Fitting Models 

In this section we describe the models used in our anal- 
ysis. We fit all the data (3 - 25 kcV) with an absorbed 
disk-corona model employing two different disk models, 
one for the whole dataset and one applied only to a sub- 
set. Details are as follow: 

3.1. Absorption model 

Wc model the absorption by the interstellar (as 
well as circumstcllar) medium by using the Tuebingen- 
Boulder interstellar medium absorption model ('tbvarab' 
in XSPEC11). The constraint on the elemental abundance 
in the direction of GRS 1915+105 is adopted from the 
recent study by Ueda et al. (2009). By using data 
from Chandra and RXTE, they make a robust constraint 
of the interstellar (plus circumstcllar) absorption compo- 
nents towards GRS 1915+105. This is, in a sense, an 
improvement in the study of GRS 1915+105 since the ab- 
sorption model significantly affects the energy range of 
the emission from a black hole accretion disk. The hy- 
drogen column density is fixed at 2.78 x 10 22 cm~ 2 with 
metal abundances larger than solar assumed. This value is 
smaller than the commonly adopted value for this source 
in which solar abundances are assumed, 4.7 x 10 22 cur 2 
(Chaty et al. 1996). 



3.2. Disk models 

We fit the disk component of the whole dataset with 
a multicolor disk (MCD) model (Mitsuda et al. 1984; 
Makishima et al. 1986), a mathematical approxima- 
tion model of the standard accretion disk (Shakura and 
Sunyacv 1973). 

Given that some of our data have a luminosity close 
to or even above the Eddington luminosity, we also fit a 
subset of the whole dataset with another disk model, the 
extended disk blackbody (extended DBB) model (the so- 
called p-free disk model; Mineshige et al. 1994). This sub- 
set contains data with disk-dominated spectra (see section 
4.1.2 for details). 

The standard theory of accretion disks gives an effective 
temperature profile which is proportional to r -3 / 4 , where 
r is the disk radius. As the luminosity approaches the 
Eddington luminosity, the standard accretion disk model 
breaks down and should be replaced by the supercriti- 
cal accretion model (e.g. the slim disk model). In the 
supercritical accretion regime, the temperature profile is 
expected to become flatter at small radii; this is caused 
by the advective motion of the trapped radiation. In the 
extended DBB model, the effective temperature profile is 
assumed to be T c ff oc r~ p , where p is now a fitting pa- 
rameter. In principle, by using the extended DBB model, 
we can distinguish the slim disk from that of the stan- 
dard disk by spectral fitting (e.g. Vierdayanti et al. 2006; 
Okajima et al. 2006; Tsunoda et al. 2006). Furthermore, 




Fig. 2. The spectral fits and residuals of some representative data: dataset 16 (left) and datasct 19 (right). The total model is 
shown in green, while the disk and Comptonization components are shown in magenta and blue, respectively. 



mass loss due to radiation pressure-driven outflow is also 
known to modify the temperature profile (Shakura and 
Sunyacv f973; Poutanen et al. 2007; Takcuchi ct al. 
2009). For a mass accretion rate of M oc r s , where s is 
a positive parameter, the effective temperature profile fol- 
lows T cff ocr~( 3 - s >/ 4 . 

3.3. Thermal Comptonization model 

We add a thermal Comptonization model to the disk 
models mentioned above ('nthcomp' in XSPEC11; Zycki 
et al. 1999) to fit the high energy part of the spectrum 
(> lOkcV), that is, the coronal component. Most of our 
data prefer low electron temperature values (< 10 keV) 
and nthcomp is one of the better Comptonization mod- 
els for such cases. It is also a consistent model for our 
present study since the seed photons can have a disk spec- 
trum (compare, for example, with 'comptt' in XSPEC11 in 
which the seed photons are assumed to have a Wien spec- 
trum). The temperature of the seed photons of the corona 
is linked to the disk temperature, by assuming that the 
seed photons are coming from the disk. 

To make the fitting realistic, we also consider the reflec- 
tion of up-scattered photons by the disk based on the code 
by Magdziarz & Zdziarski (1995) (pexriv in XSPEC11). 
Here, however, we have to apply several assumptions. We 
assume that the disk cmissivity profile is proportional to 
r -3 , where r is the disk radius. Other important parame- 
ters for the reflection component are the reflection ampli- 
fication (= J7/27r), the ionization parameter, £ (= L/nr 2 , 
where L is the luminosity and n is density) , and the inner 
and outer radius of the reflecting area of the disk, R ln and 
-Rout , respectively. The iron-K emission line is included in 
the model. 

Throughout our study we treated these reflection pa- 
rameters as free parameters except for i? ut that is always 
fixed at 10 5 r g (r g = GM/c 2 ). We find that most of our 
data prefer low R ln (< 100r g ), inferring significant reflec- 
tion from the inner part of the disk although its error is 



large. We thus fix R in at 40r g , adopting the typical value 
from the previous study by Done et al. (2004). Low values 
of £ are also preferred by our data. Therefore, we adopt 
the value of £ as found in Ucda ct al. (2009); that is, we 
fix £ at 40. The overall fitting results, however, do not 
sensitively depend on our choice of £. Finally, we treat 
J1/27T as the only free parameter of the reflection compo- 
nent. We, however, limit the maximum value of £l/2ir to 
1. 

Several studies have reported the discovery of absorp- 
tion line features from this source (see, e.g. Kotani et 
al. 2000). We find that the fitting residuals improved 
significantly when we include an iron edge (modeled with 
'edge' in XSPEC11) to account for absorption features of 
Fe XXVI ions as done in Ueda et al. (2009). We thus 
expect that H-like absorption line of iron to be present 
in the spectra of our final sample. Therefore, we include 
an absorption line at 7.0 keV (Ka). Regarding the ab- 
sorption lines, we use a negative Gaussian model. We fix 
the line width at 40 eV but left the normalization of this 
model as a free parameter. In addition, this normalization 
is also linked to the absorption depth of the iron to make 
a more physical fitting. This is necessary because both 
the equivalent width of an absorption line and the optical 
depth of an edge are proportional to the normalization 
of the Gaussian, provided that the absorption equivalent 
width is in a linear part of the curve of growth (see Kotani 
et al. 2000). 

4. Fitting Results and Analysis 

As mentioned earlier, the energy range for fitting is 
set to be 3 - 25 keV (PCA data only) and there are 25 
datasets in our final sample (see table 1). We first fit all 
the datasets with the (MCD+nthcomp) model. See Figure 
2 for two representative samples of disk-dominated spec- 
tra (right panel) and Comptonization-dominated spectra 
(left panel). 



GRS 1915+105 at Near-Eddington Luminosity 



5 




Soft color 



Fig. 3. Color-color diagram of our sample. The numbers rep- 
resent the order of the observation time listed in table 1. Most 
data inside the blue circle have Comptonization-dominatcd 
spectra and those inside the red circle have disk-dominated 
spectra. 

A color-color diagram of our data sample is given in fig- 
ure 3. We follow Done and Gicrlinski (2003) in defining 
the soft and hard colors. That is, the soft color is defined 
as the ratio of unabsorbed fluxes in the 4 - 6.4 and 3 - 
4 keV bands, while the hard color as the ratio of unab- 
sorbed fluxes in the 9.7 - 16 and 6.4 - 9.7 keV bands. As 
mentioned in section 2, we also include two datasets which 
were previously studied by Done et al. (2004). Datasets 1 
and 2 in this work correspond to datasets 6 and 7 of Done 
et al. (2004). The upper left group in figure 3 consists 
of data with Comptonization-dominated spectra with one 
exception, dataset 25. On the other hand, the lower right 
group mostly consists of data with disk-dominated spec- 
tra with two exceptions, datasets 11 and 12. Dataset 2 
is Comptonization-dominatcd when fitted with our model 
as opposed to the previous study. 

A second disk model (the extended DBB model) was 
applied to the subset in the lower right of the color-color 
diagram (see figure 3, red circle) and dataset 25 from the 
upper left group. In total, we fit 12 datasets with the 
(extended DBB modcl+nthcomp) model. We obtain ac- 
ceptable values of reduced chi-squared, x 2 /d.o.f. where 
d.o.f. means degrees of freedom, for overall fitting (ex- 
cept dataset 2, where x 2 /d.o.f.> 1.5). In fact, in most 
cases, the x 2 /A.o.i. is rather low. We found that this 
is caused by the moderately large systematic error value 
(1.5%) included in the fitting. The results from fitting 
all the data with the (MCD+nthcomp) model are given 
in tables 2 and 3, while the results from the (extended 
DBB+nthcomp) model are given in table 4. The unfolded 
spectra of all datasets in our final sample, together with 
the corresponding components of the models, are given in 
figures 4 and 5 for the (MCD+nthcomp) model and figure 
6 for the (extended DBB+nthcomp) model. 



Jf..l. Disk Components 

4- 1.1. MCD components 

Our final sample has kT ln ranges between 0.7 - 2 
keV while the total luminosity ranges from ^0.15 - 
1.04Le, where Le is the Eddington luminosity (= 1.5 x 
10 38 (M/M Q ) erg s^ 1 ). Data with kT- m within 1.5 - 2 
keV are mostly those with disk-dominated spectra (see 
figure 7, filled stars). The inner disk temperature in 
our sample is too high for a 14M non-rotating black 
hole. A rotating black hole is, thus, often suggested to 
explain a too high disk temperature phenomena in this 
source (McClintock et al. 2006; Middleton et al. 2006; 
Zhang et al. 1997). Caution needs to be taken for 
some Comptonization-dominated data. As shown in ta- 
ble 2, some of these Comptonization-dominated data (i.e. 
datasets 16, 17, 18, 20, 21, and 23) have a very low disk 
fraction , ~ 0. In this case, the values of the disk temper- 
ature from the fitting become less reliable since the disk 
contribution is negligible. 

The disk inner radius can be derived from the nor- 
malization of the MCD model, i.e., normalization = 
((£>/10kpc) 2 /(rj n [km]) 2 ) x cosi, where D is the distance 
to the source in kpc and i is the inclination angle of the 
system. Note that r ln is the apparent inner disk radius 
without relativistic (£) and spectral hardening factor (k) 
corrections. The true inner radius can be obtained by us- 
ing the following formula: r\™° — K 2 £r?f parent (see Kubota 
et al.1998). That is, when we mention r- m , we refer to 
the r*™ e throughout this paper unless stated otherwise. 
However, since we use a disk-corona model, we have to 
take into account the effect of the thermal Comptonization 
in the spectrum. We calculate the apparent disk inner 
radius by using the formula (Al) shown in Kubota and 
Makishima (2004) in which spherical geometry is assumed 
for the thermal Comptonization emission. However, here 
we assume that half the photons in the corona are injected 
back into the accretion disk due to the large optical depth, 
as done in Ueda et al. (2009). 

By using this formula, the obtained true inner radius of 
the disk is mostly within ~ 0.5 - 7rs, where r§ is about 42 
km for a 14M Q black hole (r s = 2GM /c 2 is a Schwarzschild 
radius), assuming a non-rotating black hole. Here we have 
used k = 1.7 and £ = 0.412 to derive these true inner radii. 
We find that the data whose disk fraction is larger than 
50% of the total component have a low inner disk radius, 
< Irs ( sec figure 7, filled stars). If, however, we choose 
higher k, say n = 3 (note that in figure 10 the upper branch 
in the lower panel coincides with the line of k = 3 for a 
non-rotating black hole), r m would range within ~ 1.5 - 
22rs, that is (3/ 1.7) 2 times larger than that of ft = 1.7 
case. 

We find a strong correlation between the disk parame- 
ters and the fraction of the disk component. As shown in 
figure 7, the inner disk temperature increases as the frac- 
tion of the disk component increase. The inner disk radius 
follows the opposite trend. We attempt to interpret this 
correlation in section 5. 
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Fig. 4. Unfolded spectra and model of datasets 1 - 12. The total, disk (MCD), and thermal Comptonization components are 
in green, magenta, and blue lines, respectively. The data, shown with red crosses, coincide with the total model. 
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Table 2. MCD components of the fitting. The unit for kTi n is keV. r- ln is not calculated from the MCD model normalization (see 
text). -Fdiak an d ^jigjj are the 0.01 - 100 flux (10 -8 erg cm~ 2 s _1 ) and the photon flux (in photon s _1 cm -2 ) of the disk component, 
respectively. % Disk is the disk fraction on the spectrum in per cent. The values of -Fdiski ^disk' % Disk, and i/disk/^E which are 
less than 10~ 5 are written as 0. 



Dataset 



r in /r s F disk % Disk L disk /L E X 2 /d.oi 



1 

2 
3 
4 
5 
6 
7 
8 
9 

10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 



l-8li ;J 6 

1 06+ 01 

i - uu -0.09 



1.71 
1.74 



+0.03 
-0.03 

+0.03 
-0.02 



0.84ft" 



0.89 
2.05 



+0.07 

-0.15 

+0.02 
-0.03 



1 O7+0.04 
1-O'_0.03 



1.93 
1.92 
1.57 



+0.04 
-0.03 

+0.04 
0.03 

+0.05 
-0.04 



i-68±8:S 



1.73 



+0.03 
-0.03 



i OO+0.03 
i - ZJ -0.03 



1.26 



+0.03 
-0.03 



1.10 



+0.08 
-0.03 



1.00 



+0.12 
-0.02 



0.87 



+0.04 
-0.04 



1.84 



+0.03 
-0.03 



0.86 



+0.04 
-0.04 



0.76 



-0.09 
-0.08 



90+ 11 

u - yu -0.04 



1.27 



-0.05 



1 r 1+0.03 

i ' O1 -0.03 



0.50 ±0.07 
1.80 ±0.25 
0.57±0.08 
0.55 ±0.07 
3.45 ±0.50 
2.42 ±0.35 
0.21 ±0.03 
0.49 ±0.07 
0.49 ±0.07 
0.46 ±0.06 
1.09±0.15 
1.06 ±0.15 
0.82±0.11 
1.64 ±0.24 
1.26 ±0.18 
3.37 ±0.49 
4.12±0.59 
4.93 ±0.71 
0.48 ±0.07 
3.76 ±0.54 
4.14 ±0.60 
6.61 ±0.95 
2.72 ±0.39 
1.33 ±0.19 
0.49 ±0.07 



5.21 
0.95 
1.63 
1.68 
0.001 
0.47 
3.70 
2.57 
2.46 
2.66 
1.87 
2.48 
4.09 
2.21 
1.34 




7.43 



0.28 


1.12 
0.98 



22.8 
6.77 
7.52 
7.61 
0.009 
3.98 
14.4 
10.9 
10.1 
11.0 
9.52 
11.7 
18.6 
13.9 
8.27 




32.1 




2.81 


6.82 
5.1 



74.6 
29.0 
53.0 
53.5 
0.02 
13.7 
87.3 
62.8 
59.5 
63.0 
32.8 
36.7 
57.3 
34.9 
33.5 




82.3 



2.85 

24 
51.1 



0.68 ±0.20 
0.12 ±0.04 
0.21 ±0.06 
0.22 ±0.06 
1E-4 ± 4E-5 
0.06 ±0.02 
0.48 ±0.14 
0.33 ±0.10 
0.32 ±0.10 
0.35 ±0.10 
0.24 ±0.07 
0.32 ±0.09 
0.53 ±0.15 
0.29 ±0.08 
0.17 ±0.05 




0.97 ±0.28 



0.04 ±0.01 


0.15 ±0.04 
0.13 ±0.04 



0.22 
1.57 
0.46 
0.37 
0.52 
0.43 
0.42 
0.28 
0.51 
0.57 
0.20 
0.30 
0.32 
0.25 
0.24 
0.29 
0.29 
0.24 
0.43 
0.38 
0.37 
0.26 
0.52 
0.40 
0.40 
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Table 3. Thermal Comptonization component of the fitting. The unit for kT e is keV. -Fthc an d ^xhc are 0-01 ~~ 100 keV flux 
(10 — 8 erg cm -2 s — 1 ) and the photon flux (in photon s^ 1 cm -2 ) of the Comptonizcd component, respectively, t is the scattering 
optical depth and L is the total luminosity (see text). 



Dataset kT c V F TH c F% HC T <= ^/ 2lT 

1 3.08±g;g| 2.12±g;g| 1.77 4.83 9.00 0.69 0.76 ±0.22 

2 3.34±g:§g 3.16±g:g| 2.31 12.6 5.09 1.00 0.23 ±0.06 

3 10.6li.s7 3 - 74 io:o6 I- 44 5.35 I- 88 4E " 8 0.28 ±0.08 

4 ll.Si^oi 3.75ig;g? 1.46 5.29 1.78 2E-6 0.28 ±0.08 

5 i7.iii:?o 2 - 76 io:oi 4 - 25 23.0 2.10 °- 30 0.19 ±0.05 

6 25.0i^ 3 2.58ig;g? 3 2.94 13.4 1.78 0.41 0.19 ±0.06 

7 3.38±g:gf 1.94±g;g| 0.54 1.22 9.64 1.00 0.51 ±0.15 

8 4 -60io;24 3.27ig;o6 1-52 5.10 3.99 0.41 0.40 ±0.12 

9 6.30i£s! 3.8lig;g| 1.67 5.71 2.64 0.30 0.39 ±0.11 

10 3 -61±aio 2 - 95 io:o4 I- 56 5.00 5.31 °- 49 0.42 ±0.12 

11 3 - 38 iao5 3 -09i°;o2 3 - 84 15.0 5.19 0.21 0.41 ±0.12 

12 5.71+8:46 4 - 4 0io;o9 4 -29 17.4 2.32 0.27 0.51 ±0.15 

13 5.53±g;|| 3.98ig;g? 3.05 11.5 2.72 0.37 0.67 ±0.19 

14 25i^ 4 2.89ig;$j 4 4.12 16.4 1.49 0.12 0.47 ±0.14 

15 16.6i2:?o 2 - 86 io:o2 2.67 1°- 5 2.03 0.21 0.29±0.08 

16 2 5i°:87 3 - 33 iao2 4 9 - 94 49 - 8 1-18 0.24 0.44 ±0.13 

17 14.4ii;p 3.15±g;g? 10.6 56.2 1.95 0.27 0.47 ±0.14 

18 H-7ia?2 3.04ig:81 9.09 53.3 2.37 0.31 0.40 ±0.12 

19 7.51±g;g§ 2.84i°;^ 1.60 4.77 3.52 0.51 1.04 ±0.30 

20 H- 2 ia68 3.03ig:8i 5.01 29.7 2.46 0.32 0.22 ±0.06 

21 15.6±i|g 3.20±g;gl 14.4 71.6 1.80 0.27 0.64±0.18 

22 10.6io:g4 2 - 9 lio:oi 9 - 66 62 -l 2.70 0- 42 0.47 ±0.13 

23 25+^0! 2.85±g:gl 3 - 47 17 - 9 L51 0.29 0.15 ±0.04 

24 14.lii:2? 2.52i°;°i 3.54 12.0 2.76 0.12 0.30 ±0.09 

25 21.2±|'|| 2.50i°;^ 0.94 2.10 2.69 4e-5 0.17 ±0.05 
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Fig. 6. Unfolded spectra and model of the disk-dominated data. The total, disk (p-free disk), and thermal Comptonization compo- 
nents are shown in green, magenta, and blue lines, respectively. The data, shown with red crosses, coincide with the total component. 
See text for details. 
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Table 4. Fit results of the extended DBB plus thermal Comptonization model. The units for kT[ n , kT c , Ldisk/^E, and L/Le are 
the same as in table 2 and 3. 



Dataset 



kT u 



P 



kT c 



% Disk i disk /L E 



L/L E X 2 /d-o.f 



1 

3 
4 
7 
8 
9 

10 
11 
12 
13 
19 
25 



1.85 
2.06 
2.11 
2.03 
1.91 
2.06 
1.94 
2.14 



-11 


03 


'-0 


.03 


+ 


.01 


'-0 


01 


+ 


.02 


-0 


.02 


+0 


.02 


'-0 


02 


+ 


.02 


-0 


02 


i n 


no 

\)Z. 


>±o 


.02 


+0 


.02 




02 


+ 


03 


'-0 


01 


+ 


15 


'-0 


.15 


, + 


,05 


-0 


01 


,+0 


14 


-0 


14 


+0 


12 


'-0 


04 



0.69 ±0.12 3.64 ±0.76 2.93 ±2.00 

0.59 ±0.02 8.09 ±1.86 2.86 ±0.77 

0.59 ±0.02 13.6 ±3.60 2.93 ±0.76 

0.74 ±0.05 8.28 ±1.86 4.43 ±0.75 

0.72 ±0.04 7.05 ±1.15 4.05 ±0.05 

0.68 ±0.03 25.0 ±0.0 4.52 ±0.09 

0.75±0.06 5.13±0.11 3.94±0.05 

0.56±0.01 4.92±0.07 3.15±0.03 



56.7 0.55 ±0.16 0.70 ±0.20 0.45 

85.0 0.48 ±0.14 0.51 ±0.15 0.92 
87.2 0.49 ±0.14 0.51 ±0.15 0.79 
54.2 0.30 ±0.09 0.39 ±0.11 0.67 

54.4 0.30 ±0.09 0.39 ±0.11 0.64 
67.9 0.40 ±0.11 0.46 ±0.13 0.78 
48.7 0.27±0.08 0.36±0.10 0.67 

83.1 1.06 ±0.30 1.13±0.33 1.03 
87.6 1.06 ±0.30 1.11 ±0.32 0.69 

81.5 0.89 ±0.25 0.96 ±0.28 0.76 
88.9 1.12 ±0.32 1.17±0.34 0.28 
76.0 0.33 ±0.10 0.37 ±0.11 0.60 



4-1-2. Extended disk blackbody components 

For consistency, when the extended DBB model is used 
to fit the disk component, we modify the temperature 
profile of the seed photons in the thermal Comptonization 
model. That is, the temperature profile is set to be T e g oc 
r -p , with p ranging from 0.5 to 0.75, instead of using the 
MCD temperature profile. 

We only fit a group of data in the lower right of color- 
color diagram (see figure 3, red circle) and dataset 25 from 
the upper left group (blue circle). That is, we only choose 
data with disk-dominated (or suspected, such as datasets 
11 and 12) spectra. We are interested in finding devi- 
ations from the standard high/soft state spectra as the 
luminosity approaches or exceeds the Eddington luminos- 
ity. In total, we fit 12 datasets with the extended DBB 
modcl+modificd thermal Comptonization model. 

The unfolded spectra and all components of the model 
are shown in figure 6. Results of the fits are summarized 
in table 4. As shown in table 4 (and figure 6), the disk 
fraction of the twelve datasets is mostly larger than 75%. 
Moreover, disk fraction in datasets 11 and 12 increases 
significantly and even exceeds the fraction of the thermal 
Comptonization component. 

The radial temperature gradient, p, is shown in column 
3 of table 4. We can see that the p- values tend to deviate 
from 0.75. In fact, nine out of twelve datasets (1, 3, 4, 9, 
11, 12, 13, 19, and 25) prefer a p- value which is less than 
0.7 with less than 20% error. Moreover, four out of the 
nine datasets (3, 4, 11, 25) prefer a p-value which is less 



than 0.6 with less than 10% error. The disk temperature 
of most datasets tends to increase compared to those of the 
MCD+thermal Comptonization model, that is kT ln ~ 1.8 
-2.1 keV. 

4-2. Thermal Comptonization Component 

Here we will only focus on two important parameters 
from the thermal Comptonization model that we use in 
our present study. These parameters are the electron tem- 
perature, kT D , and the power-law photon index, V . The 
electron temperature determines the high energy rollover. 
In addition, the low energy rollover is also considered and 
is parametrized by the seed photon temperature. In our 
present study, this seed photon temperature is equal to 
the inner disk temperature. The power-law photon in- 
dex determines the slope between low and high energy 
rollovers. 

We will first describe the results for the case where the 
disk component is modeled by MCD. With some few ex- 
ceptions, data with > 50% disk fraction prefer low electron 
temperature below 10 keV (see figure 8, filled stars). Wc, 
however, set the maximum limit at 25 keV for the electron 
temperature since, otherwise, it cannot be constrained by 
our data in some cases. It is important to note here that 
the error estimate for the electron temperature is, unlike 
other parameters, quite large. In other words, we should 
be very careful in making an interpretation based on the 
electron temperature results. The power-law photon in- 
dex varies between 2 - 5 for most data regardless of the 
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Fig. 7. fcT; n and r; n /rg obtained from the fitting. The 
dotted line in the lower panel shows fcT; n = 1 keV. The 
dotted lines in the upper panel mark r; n /rg =1 and 3. 
Comptonization-dominated data are shown as open squares 
while disk-dominated data are shown as filled stars. 



Fig. 8. fcT e and T obtained from the fitting (see text). The 
dotted line in the lower panel represents kT c = 10 keV. 



disk fraction. Despite its relatively small errors, no clear 
trend is found which left no room for further speculation 
(see figure 8). 

In figure 9, the electron scattering optical depth, r 
is calculated from T c and T by following Sunyaev and 
Titarchuk (1980): 



.25 



(T e /511kcV)[(r + 0.5) 2 -2.25] 



-1.5. (1) 



The lines show constant T = 2, 2.5, 3, 3.5, 4, and 4.5, 
from top to bottom. All data prefer an electron scattering 
optical depth between 1 and 10, mostly greater than 2. 

In the case where the disk component is modeled by ex- 
tended DBB, the electron temperatures tend to be higher 
but mostly still less than 10 keV, in agreement with those 
of the MCD fits (see tables 3 and 4) . The power- law pho- 
ton index, T, remains within the same range as those of 
the MCD fits although the values are not always similar. 
We conclude that both the MCD and the extended DBB 
model prefer low temperatures (< 20 keV) and thus opti- 
cally thick (r > 2) corona (since T values do not vary so 
much). 




(keV) 



Fig. 9. t vs. fcT e (see text). The lines are for constant value 
of r, that is T = 2, 2.5, 3, 3.5, 4, 4.5, from top to bottom. 
The open squares are data with Comptonization-dominated 
spectra, while the filled stars are those with disk-dominated 
spectra. 
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Fig. 10. Luminosity vs. disk temperature diagrams. We as- 
sumed d= 12.5 kpc and i = 70°, and M = 14M Q . The lines in 
each panel represent k = 1.7, 3.0, and 4.0 from top to bottom, 
respectively, for a non-rotating black hole (see text). Open 
squares are those with disk fraction < 50% and filled stars for 
> 50%. 

4-- 3. X-ray Hertzsprung-Russell (HR) Diagram 

One of the most general properties of accreting black 
hole systems is that they follow aloe T* s -track (see, e.g. 
Gierlihski and Done 2004) in a luminosity versus inner 
disk temperature diagram (an X-ray Hertzsprung-Russell 
(HR) diagram) during the high/soft state (i.e. when the 
disk component dominates the soft X-ray spectrum). This 
property means that as the disk temperature and luminos- 
ity change, the inner radius of the disk remains constant 
(Makishima et al. 1986; Ebisawa et al. 1991; Dotani et al. 
1997; Kubota et al. 2001). Thus, the inner disk radius is 
often related to the innermost stable circular orbit (ISCO) 
which in turn can be related to the black hole mass. 

However, as the luminosity of the system approaches 
its Eddington luminosity, it starts to depart from the L oc 
T^ s -track. Some systems, indeed, have shown that they 
start to follow a new flatter track of L oc T c 2 ff (Kubota and 
Makishima 2004), which satisfies the slim disk prediction 
for a non-rotating black hole. 

Figure 10 shows the X-ray HR diagram of our sample. 
Here we use T m instead of T e ff, but as long as we consider 
a constant spectral hardening factor, k, T- m oc T c ff ■ The y- 
axis of the upper panel shows the total luminosity (L), ob- 



tained by summing the luminosity of the disk component 
(^disk) and thermal Comptonization component (Lthc), 
normalized by the Eddington luminosity. The lower panel 
shows the disk luminosity normalized by the Eddington 
luminosity. The disk luminosity is calculated from the un- 
absorbed flux of MCD component between 0.01 - 100 keV 
assuming a disk geometry, a distance of 12.5 kpc, and an 
inclination angle of 70°. We assume an isotropic emission 
when wc calculate the thermal Comptonization luminos- 
ity (also shown in table 3). The filled stars show data with 
> 50% disk fraction in their spectra, while open squares 
represent the data with lower disk fraction (< 50%). Data 
with Ldisk < 0.02Le are not shown in the lower panel for 
clarity. 

It is not always clear which luminosity (L or £disk) 
is better for this kind of plot. In the present study 
we plot both luminosities for completeness. In disk- 
dominated data, the disk fraction is moderately higher 
than that of the thermal Comptonization component and 
thus both luminosities give similar results. By contrast, 
for Comptonization-dominated data, their position in the 
X-ray HR diagram depends on the choice of luminosity, 
L or Ldisk (compare upper and lower panels of figure 10). 
As shown in the upper panel, these data have compa- 
rable total luminosities to those of disk-dominated data. 
However, their disk fraction is mostly low and thus their 
position is shifted to a lower position in the lower panel of 
figure 10. Data with a very low disk fraction (~ 0) do not 
follow the Ld oc X^-track in the lower panel. Moreover, 
the disk temperature value obtained from fitting becomes 
less reliable since the disk contribution is almost negligible 
(see section 4.1.1). 

From the lower panel of figure 10, we can see that 
there are two groups of data following two different Ld oc 
7]t-tracks. The upper branch seems to extend to the 
lower disk luminosity data (open squares) but these data 
are those whose spectra are dominated by the thermal 
Comptonization. All spectra of the data in the lower 
branch are disk-dominated and their disk luminosity is be- 
low 0.5-Le- The lower branch may correspond to the high- 
est end of standard accretion disk branch in the high/soft 
state or slim-disk state (see below). 

In figure 10, we also plot some lines from equation 
[9] of Makishima et al. (2000), assuming Le = 1.5 x 
10 38 (M/M Q ) erg s -1 . It is assumed that the inner edge 
of the disk reaches the ISCO, r- ln = 3r$ for a non-rotating 
black hole and r; n = 0.5rs for a maximally rotating black 
hole. Given that wc know the black hole mass from dy- 
namical methods, M = 14M Q , we can plot the lines for 
different values of k, by assuming a value of black hole 
spin. In figure 10, we assume a black hole spin of (non- 
rotating black hole). These lines of constant k coincide 
with L oc l£ -tracks. The upper and lower branches, in 
the lower panel of figure 10, coincide with n ~ 3 and 4 for 
a non-rotating black hole, respectively. 

Figure 11 shows the X-ray HR diagram for both MCD 
(blue points) and extended DBB fits (red points). The 
blue open squares and filled stars are Comptonization- 
dominated and disk-dominated data, respectively. Red 
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Fig. 11. X-ray HR diagrams for both MCD and extended 
DBB fits. Blue colors are from MCD fits: open squares for 
Comptonization-dominated data and filled stars for disk-dom- 
inated data. Red colors are from extended DBB fits: p ~ 0.75 
shown as open circles, while p < 0.7 as filled circles. Open di- 
amonds show data from full-relativistic slim disk simulations 
where ft is fixed at 1.7. The lines are L oc Tj* -lines which cross 
each simulation data point (see text). 

filled circles are those datasets with p < 0.7 from extended 
DBB fits, while red open circles are those with p ~ 0.75. 
We also plot the fitting results of two datasets from numer- 
ical simulations of the slim disk model around a rotating 
black hole (Sadowski 2009) in blue open diamonds. The 
higher temperature data corresponds to a black hole spin 
of 0.99 while the lower temperature to a black hole spin of 
0.52. Relativistic effects are also included when calculat- 
ing the spectrum and the spectral hardening factor is 1.7 
(Sadowski et al. 2009). We draw lines that cross each nu- 
merical simulation point (blue open diamond) assuming 
L oc in orange and green. 

It turns out that the dependence of L on Ti n becomes 
less unclear because of the limited range of T m . Caution 
needs to be taken, however, since the luminosity calcula- 
tion from the observed flux for the extended DBB model is 
less reliable than that of the MCD model because the disk 
state is expected to become the standard one and hence 
p should change to 0.75 at outer radii. Rather, we can 
directly compare our MCD fit results (blue points) with 
the same plot obtained from other BHBs (e.g. Figure 5 
of Kubota and Makishima 2004) , where the L oc -track 



is suggested at the highest temperature end based on the 
MCD model. 

5. Discussion 

GRS 1915+105 is known to be the only source in 
our Galaxy that persistently shines at above 10% of its 
Eddington luminosity. It raises a doubt on whether we 
ever observed the canonical high/soft state even when the 
disk fraction is large. We also found a puzzling behavior 
from our present study. In this section we will discuss our 
interpretation as well as explore the available scenarios for 
this puzzling behavior at near-Eddington luminosity. 

5.1. Accretion Disk States 

As shown in figure 7, there is a clear dependence of the 
MCD fitting components, i.e. rj n and kT m on the disk 
fraction. As the disk fraction increases from to 50%, 
the innermost temperature of the disk increases as the 
radius decreases. It is also shown that n n remains almost 
constant for those datasets with disk fraction larger than 
50%. 

We suggest that the Comptonization-dominated spectra 
most probably correspond to the (strong) very high state 
(VHS) found in other BHBs, because of the similarity in: 
(1) the parameters of Comptonization (kT c ~10-20 keV 
and T e ~ 2-3), (2) the high L/Le value, (3) the (appar- 
ently) large r- ln (or high L) for the given T in , and (4) the 
presence of strong (low frequency) QPOs (e.g. Morgan et 
al. 1997; Rodriguez et al. 2004). These are exactly found 
in the strong VHS of XTE J1550-564, although the real 
innermost radius may be smaller if the disk-corona cou- 
pling is taken into account (Done and Kubota 2006). This 
VHS occurs when the mass accretion rate becomes even 
higher than that of the high/soft state. Theoretically, 
some outflow material blown off the disk due to the in- 
crease of radiation pressure caused by the increase of mass 
accretion rate is also predicted in this state. It is very 
unlikely that the Comptonization-dominated spectra cor- 
respond to the canonical low/hard state of BHBs, which 
is normally observed at luminosities much lower than 10% 
of the Eddington luminosity, unlike our case. In addition, 
the electron temperature of the Comptonizing corona in 
the low/hard state (e.g. Makishima et al. 2008) is much 
higher (~100 keV) than those found here. 

We propose that the following three branches (states) 
have been observed in our data sample: 

1. Comptonization-dominated spectra (when the disk 
flux < 30% of the total flux), which correspond to 
the strong VHS in other BHBs. 

2. Disk-dominated spectra: 

(a) Lower luminosity branch: the high/soft state, 
or probably the slim disk state, around a ro- 
tating black hole (see next subsection). 

(b) Higher luminosity branch: probably a trun- 
cated disk with low temperature (< 10 keV) 
and optically thick Comptonization. This 
would correspond to the new state that appears 
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Fig. 12. Cartoon of the three branches from the upper panel 
of figure 10. 

at luminosity very close to the Eddington lumi- 
nosity. It may be understood as the extension 
of the high electron temperature VHS (state 1 
above) when luminosity becomes higher. 

State (1) and (2a) have already been observed in other 
BHBs, but state (2b) looks unique in GRS 1915+105, 
probably due to its high L/Le- The disk-dominated spec- 
tra (state 2) will be discussed in more detail in the fol- 
lowing subsection. The cartoon of these three branches 
in the X-ray HR diagram is shown in figure 12, which 
corresponds to the upper panel of figure 10. 

5.2. On Disk- dominated Spectra and the Black Hole Spin 

The fact that the L oc T^-track does not appear in the 
disk-dominated branches even when the luminosity ap- 
proaches and even exceeds the Eddington limit contra- 
dicts the theoretical predictions of slim disk model for a 
non-rotating black hole. Moreover, as shown in the lower 
panel of figure 10, a high value of the spectral hardening 
factor, k ~ 4 is required by the data in the lower branch 
assuming a non-rotating black hole. The higher value of 
k compared to the standard value commonly used for the 
standard case k = 1.7 may be related to the increase in 
the mass accretion rate. 

Beloborodov (1998) has calculated the disk structure 
of a rapidly rotating black hole with a super-Eddington 
accretion rate. The calculations show that when the vis- 
cosity is high and the mass accretion rate approaches the 
critical value, a strong deviation from equilibrium occurs 
in the inner region of the disk. In this region, the plasma 
density is low and thus emission capability is also low. 
The produced radiation in the disk is thus not fully ther- 
malized which results in higher disk temperature which 
further can be related to higher value of k. The ques- 



tion remains, however, since the luminosity of the lower 
branch, whose K is higher, is lower than that of the upper 
branch. 

We may have underestimated the luminosity of the 
lower branch. For example, if some part of the disk is ob- 
scured by highly ionized outflow materials which scatter 
disk radiation. In the case when the outflow has a small 
opening angle, we should have overestimated the luminos- 
ity due to high inclination angle of GRS 1915+105 and not 
the other way around. If, however, we have a moderately 
geometrically thick disk (slim disk) together with an out- 
flow, there may be some interaction between the disk and 
the outflow and complicate the radiation track. The ma- 
terial outflow from the disk is a delicate process which far 
from being understood. We need a numerical calculation 
of disk-corona system with an outflow for further analy- 
sis. Kawashima et al. (2009) performed two dimensional 
radiation hydrodynamic simulations of supercritical accre- 
tion flows with outflows. They found a possible new state 
above the slim disk state in which both disk and outflow 
components become important. The slim disk and out- 
flow model seems to be applicable in a particular case of 
ULXs (Gladstone et al. 2009). 

High inner disk temperatures found in our study may 
be explained by a slim disk around a non-rotating black 
hole (e.g. Watarai et al. 2000), a standard disk around 
a rotating black hole (e.g. Zhang et al. 1997), or a slim 
disk around a rotating black hole (e.g. Sadowski 2009). 
As mentioned above, the non existence of L oc -track 
and the too high k value contradict the theoretical pre- 
diction of the slim disk around a non-rotating black hole. 
While the rotating black hole is a good explanation for 
the high inner disk temperatures (since the inner disk ra- 
dius becomes closer to the black hole) , the luminosities of 
our sample are too high to keep the disk in the standard 
accretion mode. 

On the other hand, a clear deviation from the standard 
accretion disk spectra during the high/soft state as lumi- 
nosity approaches the Eddington luminosity is in a good 
agreement with slim disk prediction. Slim disk around 
a rotating black hole is one plausible idea to resolve the 
above issues since the L oc -track is not necessarily held 
for the following reason. 

In the case of a disk around a non-rotating black hole, 
there is a gap between the radius of the ISCO (i.e., 3rg) 
and that of the event horizon (Irs). As the mass accretion 
rate increases, the flow density inside the ISCO increases 
and eventually fills the gap, emitting significant radiation. 
Hence, substantial radiation will be emitted from there 
at luminosities near the Eddington. This results in an 
apparent decrease of the inner edge of the disk and in- 
crease of the innermost disk temperature with an increase 
of the mass accretion rate (and luminosity), in the way 
that r in ocT;" 1 . This leads to the relation, L oc rf n T^ oc T? n 
(Watarai et al. 2000). This cannot happen for a maxi- 
mally rotating black hole, since the ISCO is located near 
the event horizon, regardless of the mass accretion rate. 
Hence, the relation, L oc T?, will never appear. 

Furthermore, If we assume a rapidly rotating black hole, 
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we do not need to adopt too high spectral hardening factor 
for the data in the lower branch of figure 11 (lower panel). 
Provided that the spin parameter does not change, the 
nature of the upper branch remains unclear . One pos- 
sible explanation is the change of spectral hardening fac- 
tor despite the unknown mechanism. Another possible 
explanation is the effect of the strong Comptonization. 
The inner disk is truncated and turns into Comptonizing 
corona, which is suggested as a picture of the VHS of nor- 
mal BHBs, and thus we observed lower inner disk tem- 
perature than expected from high/soft state at the same 
luminosity (Done and Kubota 2006). 

The fact that the upper branch in our study has never 
been observed in other Galactic microquasar, such as XTE 
J1550-564 and GRO J1655-40 (e.g. Kubota et al. 2001; 
2004), may have something to do with the spin parameter. 
Both XTE J1550-564 and GRO J1655-40 are suggested 
to have a moderate spin parameter < 0.8 (e.g. Gierlihski 
etal. 2001; Davis ct al. 2006; Shafee et al. 2006). Caution 
needs to be taken since there are several methods in de- 
termining spin parameter which can lead to contradicting 
results. In the case of GRO J1655-40, for example, Reis 
et al. (2009) obtained a lower limit for the spin of 0.9 
by using the reflection component in the spectrum. At 
present, there is no consensus on the best methods to de- 
termine the black hole spin parameter. Despite the fact 
that our present study does not focus on the determina- 
tion of the spin parameter, our results suggest a rapidly 
rotating black hole for GRS 1915+105. 

McClintock et al. (2006) find a spin parameter > 0.98 
for GRS 1915+105 from their analysis. They also find 
that the spin decreases as the luminosity increases. They, 
however, choose only data with < 30%Le to deter- 
mine the spin parameter. The spin parameters obtained 
from data with La > 40%Le, are indeed in agreement with 
the moderate spin parameter obtained by Middleton et 
al.(2006), spin parameter - 0.7. Middleton et al. (2006), 
however, due to their selection criteria, have only a re- 
stricted data sample which contains only high luminosity 
data, La > 40%-Le- Comparing the disk luminosities and 
temperatures obtained from both Middleton et al.(2006) 
and McClintock ct al. (2006), data of Middleton et al. 
(2006) seem to correspond to our higher branch, while 
those of McClintock et al. (2006) with L d < 30%L E cor- 
respond to the lower branch in our X-ray HR diagram. 

6. Conclusions 

We analyzed the spectra of GRS 1915+105 during its 
quasi-steady states to study the supercritical accretion 
processes. We found new interesting features as follows: 

1. There is a strong dependence of the inner disk tem- 
perature and radius on the fraction of the disk com- 
ponent: the disk temperature increases as the disk 
fraction increases, while the opposite trend is found 
for the inner disk radius. The electron temperature 
tends to be higher when the disk fraction is small 
while the spectral slope remains more or less con- 



stant. 

2. We also fit the data with the extended DBB model 
to see any deviation from the standard spectral 
shape. Despite the fact that the p-value deviates 
from the standard value p = 0.75 for data with 
L > 0.3Le, the L oc T^-track predicted by the slim 
disk theory for a non-rotating black hole does not 
appear. A rapid spin of the black hole in GRS 
1915+105 can explain these features as well as a 
high value of spectral hardening factor (~ 4) that 
would be required for a non-rotating black hole. 

3. We observed three spectral branches (states) in 
our data sample: one Comptonization-dominated 
branch and two disk- dominated branches. We sug- 
gest that Comptonization-dominated spectra cor- 
respond to the canonical VHS observed in nor- 
mal BHBs. The lower luminosity branch of disk- 
dominated spectra seems to correspond to the high- 
est end of the high/soft state or the slim disk state 
around a rotating black hole. The higher luminos- 
ity branch of disk-dominated spectra may be inter- 
preted as a truncated disk with low temperature 
and optically thick Comptonization, unique to GRS 
1915+105 at near Eddington luminosity. 
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